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Abstract 



Wei make a detailed measurement of the metal abun- 



ducing the bulk of the cluster optical light (Arnaud et al 
1994 |Rcnzini ct al. 1993|). 



dance profiles and metal abundance ratios of the inner 
core of M87/ Virgo observed by XMM-Newton during the 
PV phase. We use multi temperature models for the in- 
ner regions and we compare the plasma codes APEC and 
MEKAL. We confirm the strong heavy elements gradient 
previously found by ASCA and BeppoSAX, but also find a 
significant increase in light elements, in particular O. This 
fact to gether with the constant O/Fc ratio in the inner 9 



While the origin of the metals observed in the ICM 
is clear (they are produced by supernovae, in particular 
SNII and SNIa), less clear is the transfer mechanism of 
these metals to the ICM. The main mechanisms that have 
been proposed for the metal enrichment in clusters are: en- 
richment of gas during the formation of the proto-cluster 



(Kauffman & Chariot 1998); ram pressure stripping of 
meta l enriched gas from cluster galaxies (Gunn & Gott 



arcmir j indicates an enhancement of contribution in the 
core o f the cluster not only by SNIa but also by SNII. 



1972| , [Toniazzo fc Schindlcr 2001|) ; stellar winds AGN- or 



SN-induccd in Earl y-type galaxies flMatteucci fc Vcttolani 



1988| , [Renzmi 1997j ) 
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Another controversial subject is the relative contri- 
bution to the metal enrichment by the two types of su- 



pcrnovac, SNIa and SNII. Mushotzky et al. (1996) and 



Mushotzky fc Lowenstein (1997)| showed a dom inance of 



& Arimoto 1997 


, Fukazawa et al. 1998, 


Finoguenov et al 
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Dupkc & White 2000a), still indicating a predomi- 



The study of metallicity in clusters of galaxies dates back 
to 1976, with the detection of the 7 keV blend line emission 
due to highly ionized iron as a strong feature in the X-ray 



spectra of the Perseus cluster (Mitchell et al. 1976). This 
was a most significant observational discovery concern- 
ing X-ray clusters, following their identification as X-ray 
sources, because it established without doubt that the pri- 
mary emission mechanism was thermal and that the hot 
intra-cluster gas contained a significant portion of pro- 
cessed gas L jvirich_had at some point been ejected from 
stars ( Barazin 1988 ). 

The X-ray emitting hot intra-cluster medium (ICM) 
of clusters of galaxies is now known to contain a large 
amount of metals: for rich clusters between red-shift 0.3- 



nance of SNII enrichment at large radii in clusters, do not 
exclude that as much as 50% of the iron in clusters come 
from SNIa ejecta in the inner part of clusters. 

Spatially resolved abundance measurement in galaxy 
clusters are of great importance because they can be used 
to measure the precise amounts of metals in the ICM 
and to constrain the origin of metals both spatially and 
in terms of the different contributions of the two differ- 
ent type of SNe (since SNIa products are iron enriched, 
while SNII products are rich in a elements such as O, Ne, 
Mg and Si) as a function of the position in the cluster. 
The first two satellites able to perform spatially resolved 
spectroscopy, ASCA and BeppoSAX have revealed abun- 



dance gradients in cD clusters (Dupke & White 2000b, Dc 



0.4 anc. the present day the observed metallicity is about 
1/3 the solar value (Mushotzky & Lowenstein 1997, A11ce| 



Grandi fc Molendi 2001), in particular M87/Virgo (Mat 



fc Fabian 1998, Fukazawa et al. 1998, Delia Ceca et al 
2000] , Ettori ct al. 2001), suggesting that a significant frac 
tion of the ICM has been processed into stars already at 
intermediate red-shifts. 

The total mass of iron is directly proportional to the opti- 
cal luminosity from elliptical and lenticular galaxies, this 
fact strongly arguing for the metals now in the ICM hav- 
ing been produced by the massive stars of the same stellar 
population to which belong the low mass stars now pro- 



sumoto et al. 199(f " Guainazzi fc Molendi 2000|) , and v aria- 



tions in Si/Fc within a cluster (Finoguenov et al 200C) and 
among clusters ( Fukazawa ct al. 1998j ). Given the diver- 
sity of metals produced by SNIa and SNII, the variations 
in Si/Fe suggest that the metals in the ICM have been 
produced by a mix of the two types of SNe. 
Moreover with the aid of observational data on many ele- 
ments, X-ray observation can discriminate between com- 
peting theoretical models for SNIa also because the central 
region of cD clusters, where the SNIa products could dom- 
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inatc in the abundance pattern, can be considered as the 
archives of an extremely large number of SNIa explosions 



( pupkc fe White 2000tj , |Finogucnov ct al 200l[ ). 

The broad bandpass large X-ray sensitivity and good 
angular resolution of XMM-Newton enables the measure- 
ment of abundances and abundance gradients of elements 
that have not been accurately studied before. In particu- 
lar the O abundance (and O/Fe ratio) is particularly im- 
portant for estimation of the SNII contribution since O is 
expected to be produced mostly by this kind of supernova. 
Limited sensitivities of pre-XMM-Newton instruments for 
O VIII lines prevent this best measure to separate the enrich- 
ment by the two types of supernovae. 

M87 is the cD galaxy of the nearest cluster and its high flux 
and close location makes it one of the most ideal target for a de- 
tailed study of the ICM abundance down to the scale of the kpc. 
Throughout this paper, we assume Ho = 50kms _1 Mpc~ 1 , 
qo — 0.5 and at the distance of M87 1' corresponds to 5 kpc. 



2. Observation and Data Preparation 



M87/ Virgo was observed with XMM-Newton (Jansen et al 



200 ]| ) during the PV phase with the MOS detector in Full 
Frame Mode for an effective exposure time of about 39 ks. De- 



tails on the observation have bee n published in Bohringer et al 



(2001) and Belsole et al. (2001). We have obtained calibrated 
event files for the MOS1 and MOS2 cameras with SASv5.0. 
Data were manually screened to remove any remaining bright 
pixels or hot column. Periods in which the background is in- 
creased by soft proton flares have been excluded using an inten- 
sity filter: we rejected all events accumulated when the count 
rates exceeds 15 cts/lOOs in the [10 - 12] keV band for the two 
MOS cameras. 

We have accumulated spectra in 10 concentric annular re- 
gions centered on the emission peak extending our analysis out 
to 14 arcmin from the emission peak, thus exploiting the entire 
XMM field of view. We have removed point sources and the 
s ubstructures which are clearly visible from the X-ray image 
(Belsole et al. 2001) except in the innermost region, where we 
have kept the nucleus and knot A, because on angular scales 
so small it is not possible to exclude completely their emis- 
sion. We prefer to fit the spectrum of this region with a model 
which includes a power law component to fit the two point like 
sources. We include only one power la w component due to th e 
similarity of the two sources spectra (Bohringer et al. 2001). 
The bounding radii are 0'-0.5', 0.5'-l', l'-2', 2'-3', 3'-4', 4'-5', 
5'-7', 7'-9', 9'-ll' and ll'-14'. The analysis of th e 4 central 



regions within 3 a rcmin was already discussed in Molendi & 



Gastaldello (2001) 



Spectra have been accumulated for MOS1 and MOS2 inde- 
pendently. The Lockman Hole observations have been used for 
the background. Background spectra have been accumulated 
from the same detector regions as the source spectra. 

The vignetting correction has been applied to the spec- 
tra rather than to the effective area, as i s customary in the 
analysis of EPIC data (Arnaud et al. 2001). Spectral fits were 



performed in the 0.5-4.0 keV band. Data below 0.5 keV were 
excluded to avoid residual calibration problems in the MOS 
response matrices at soft energies. Data above 4 keV were ex- 
cluded because of substantial contamination of the spectra by 



hotter gas emitting further out in the cluster, on the same line 
of sight. 

As discussed in 



Molendi (2001) 



there are cross-calibration 
uncertainties between the spectral response of the two EPIC 
instruments, MOS and PN. In particular for what concern 
the soft energy band (0.5-1.0 keV) fitting six extra-galactic 
spectra for which no excess absorption is expected, MOS re- 
covered the Nh galactic value, while PN gives smaller Nh by 
1 - 2 x 10 20 cm" 2 . Thus we think that at the moment the MOS 
results are more reliable than the PN ones in this energy band, 
which is crucial for the O abundance measure. For this reason 
and for the better spectral resolutions of MOS, which is again 
important in deriving the O abundance, we limit our analysis 
to MOS data. 



3. Spectral modeling and plasma codes 

The X-ray emission in cluster of galaxies originates from the 
hot gas permeating the cluster potential well. The continuum 
emission is dominated by thermal bremsstrahlung, which is 
proportional to the square of the gas density times the cooling 
function. From the shape and the normalization of the spec- 
trum we derive the gas temperature and density. In addition 
the X-ray spectra of clusters of galaxies are rich in emission 
lines due to K-shell transitions from O, Ne, Mg, Si, S, Ar and 
Ca and K- and L-shell transitions from Fe and Ni, from which 
we can measure the relative abundance of a given element. A 
measure of the equivalent width of a spectral line is a direct 
measure of the relative abundance of a given element. This 
comes from the fact that both the continuum and line emission 
are two body processes with the continuum emissivity propor- 
tional to the electron density times the proton density and the 
line emissivity proportional to the electron density times the 
density of a given element. From the definition of equivalent 
width it is easily derived that this quantity is proportional to 
the ratio between the ion and proton densities. 

In Figure |l| we show the data of the 3'-4' bin together 
with the best fitting model calculated using the MEKAL code. 
The model has been plotted nine times, each time all element 
abundances, except one, are set to zero. In this way the con- 
tribution of the various elements to the observed lines and line 
blends become apparent. In the energy band (0.5-4 keV) we 
have adopted for the spectral fitting, the abundance measure- 
ments based on K-lines for all the elements except for Fe and 
Ni, for which the measure is based on L-lines. The K-lines of O, 
Si, S, Ar and Ca are well isolated from other emission features 
and clearly separated from the continuum emission, which are 
the requirements for a robust measure of the equivalent width 
of the lines and consequently of the abundances of these ele- 
ments. The Fe-L lines are known to be problematic, because 
the atomic physics involved is m ore complicated than K-shell 
transitions (Liedahl et al. 1995, but from the very good sig- 
nal of XMM spectra and from the experience of ASCA data 



(Mushotzky et al. 1996, Hwang et al. 1997, Fukazawa et 



al 



19981) we can conclude that the Fe-L determination is reliable. 



Some of the stronger Fe-L lines due to Fe XXII and Fe XXIV 
are close to the K-lines of Ne and Mg, respectively and blend- 
ing can lead to errors in the Ne abundance and, t o a smaller 
extent, to the Mg abundance ( Liedahl et al. 1995 , Mushotzky 
|et al. 1996[ ). Also the Ni measure is difficult due to the possible 
confusion of its L-lines with the continuum and Fe-L blend. 
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We therefore apply to the central regions (inside 3 arcmin) 
three different spectral models. 

A two te mperature model (vmekal + vm ekal in XSPEC and 
model II in ( Molendi fe Gastaldello 2001) using the p lasma 
code MEKAL( [Mewe et al. 1985| , [Licdahl et al. 1995 ). The 
metal abundance of each element of the second thermal com- 
ponent is bound to be equal to the same param eter of the 
first thermal component. This model is used (e.g. Makishima 



et al. (2001) and refs. therein) as an alternative to cooling-flow 



models in fitting the central regions of galaxy clusters. 

A "fake mult i-phase" model (vmekal + vm cflow in XSPEC 



Energy (keV) 




and model III in 



Molendi fc Gastaldello (2001) 



recen t papers (Molendi & Pizzolato 2001 , Molendi & Gastaldello 



As indicated in 



Energy (keV) 



2001) this model is used to describe a scenario different from a 
multi-phase gas, for which it was written for: the gas is all at 
one temperature and the multi-phase appearance of the spec- 
trum comes from projection of emission from many different 
physical radii. A more correct description will be given by a 
real deprojection of the spectrum. 

The third model is the analogue of the vmekal two temper - 
ature model using the plasma code APEC (smith et al. 2001). 

We can't adopt an APEC analogue of the fake multi-phase 
model because the cooling flow model calculating its emission 
using APEC is still under development. Given the su bstantial 
agreement between 2T and fake multi-phase model ( Molendi 



merits, obtained by setting all element abundances to zero 
except the one of interest, calculated using the MEKAL 
code. 

Top Panel: the lines of O, Fe, Ne, Ni, Mg and Si used 
for the analysis in the 0.5-2 keV energy range. 
Bottom Panel; the lines of Si, S, Ar and Ca used for 
the analysis in the 2-4 keV energy range. 



The elemental abundances are expressed by the relative 
values to the solar abundances from |Grevesse_ fe Sauval (1998)| 



Figure\l. Data of the 'J- 4' bin and lines of the various ele- fe Gastaldello 200l|), we can regard the 2T APEC results as 



indicative also for a fake multi-phase model. 

For the spectrum accumulated in the innermost region we 
included also a power law component to model the emission of 
the nucleus and of knot A. 

For what concern the measure of metal abundances, it is 
crucial to use multi temperature models in the inner bins, 
where the evidence of multi temperature gas is stronger. In- 
fact single temperature models substantially underestimate the 
metal abundance in multi temperature spectra, as already ob- 
served in AS CA spectra of galaxies and galaxy groups (Fe bias: 



Buote 2000). In M87 the drastic abundance dec rement found 



: Grcvcssc 



that updates the commonly used values from Anders ( 
(1989)[| 

All spectral fitting has been performed using version 11.0.1 
of the XSPEC package. 

All models discussed below include a multiplicative com- 
ponent to account for the galactic absorption on the line of 
sight of M87. The column density is always fixed at a value 
of 1.8 x 10 20 cm' 2 , which is derived from 21cm measurements 



( Lieu et al. 1996 ). Leaving Nh to freely vary does not improve 
the fit and does not affect the measure of the oxygen abun- 
dance, which could have been the more sensitive to the pres- 
ence of excess absorption. The Nh value obtained is consistent 
within the errors with the 21cm value. 



The temperature profile for M87 (Bohringer et al. 2001 
shows a small gradient for radii larger than ~ 2 arcmin and 
a rapid decrease for smaller r adii. Moreover, as pointed out 
in Molendi & Pizzolato (2001) all spectra at radii larger than 
2 arcmin are characterized by being substantially isothermal 
(although the spectra of the regions between 2 and 7 arcmin are 
multi temperature spectra with a narrow temperature range 
rather than single temperature spectra), while at radii smaller 
than 2 arcmin we need models which can reproduce the broad 
temperature distribution of the inner regions. 



in the previous analysis (Bohringer et al. 2001) follows from 
the application of the oversim plified single temperature model 
(Molendi & Gastaldello 2001), as shown in Figure H 

For the outer regions (from 3 arcmin outwards) we apply 
single temperature models: vmekal using the MEKAL code and 
vapec using the APEC code. 

We use the two different plasma codes with the aim of 
cross-checking their results. As pointed out by the authors of 
the new code APEC, cross-checking is very important, since 
each plasma emission code requires choosing from a large over- 
lapping but incomplete set of atomic data and the results ob- 
tained by using independent models allows critical comparison 
and evaluation of errors in the code and in the atomic database. 



4. Results 
4.1. Abundance profiles 

The abundance measurements obtained using the two differ- 
ent plasma codes agree for what concerns Fe, Ar and Ca; they 
are somewhat different for what concerns O, S and Ni and in 
complete disagreement for Mg and Ne. The temperature pro- 
file obtained with the two codes are somewhat different in the 
inner regions (where we compare the range of temperature ob- 
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tained by multi-temperature models), while there is substan- 
tially agreement in the outer isothermal bins. 

The models using the APEC code give a systematically 
worse description than the ones using MEKAL code; we show 
some examples in Fig. ^| In the external bins the differences in 
the fit between the two codes are due to APEC over-prediction 
of the flux of Fe-L lines from high ionization states, consider- 
ing the fact that the temperature obtained by the two codes 
are nearly coincident. For the inner regions the differences be- 
tween the two codes are further complicated by the different 
temperature range they find for the best fit. In general, where 
the temperature structure is very similar, as in the innermost 
bin, the difference is as in the outer bins in the high energy 
part of the Fe-L blend, while where the temperature structure 
is different, as in the l'-2' shown in Figure |^, the differences 
between the two codes are primary due to different estimates 
of the flux of He-like Si-K line. 

We therefore choose as our best abundance profiles those 
obtained with a 2T vmekal fit for the central regions and with 
a IT vmekal fit for the outer regions. We report the abundance 
profiles obtained in this way in Fig ^ for O, Si, Fe, in Fig ^| for 
Mg, Ar, S and in Fig ^| for Ne, Ca, Ni. Abundance gradients 
are clearly evident for Fe, Si, S, Ar and Ca; a statistically 
significant enhancement is evident in the central regions for O, 
Mg and Ni, while only Ne is substantially flat. 
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4.2. Abundance ratios 

From the abundance measurements we obtain the abundance 
ratios between all the elements relative to Fe, normalized to the 
solar value. As an example we show the two more significant 
abundance ratios, O/Fe and S/Fe, in Fig. (?], together with the 
abundance ratios obta ined by models of supernovae taken by 
Nomoto et al. (1997) and rescaled to the solar abundances 



reported in Grevesse & Sauval (1998) 



An important results is the Constance in the inner 9 ar- 
cmin of the ratio O/Fe (fitting the data with a constant gives 
X 2 = 6.9 for 7 d.of. and adding a linear component does not 

MOS data 
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Figure 2. Abundance profiles for O, Si and Fe in the inner 
3 arcmin. Full circles indicate the measurements with a 
single temperature model, while empty triangles those with 
a two temperature model. Solar units are the old ones by 
Anderi & Grevess (1989). 



Figure 3. Top Panel: data for the bin 1' -0 and con- 
tributions to x 2 for a 2T model with MEKAL code, the 
X 2 /d.o.f. of the fit is 381/217, and for a 2T model with 
APEC code, the X 2 /d.o.f. of the fit is 474/218. Contribu- 
tions to x 2 are on the same scale for direct comparison. 
Middle Panel; data for the bin 3 -4' and contributions 
to x 2 for a IT model with MEKAL code, the x 2 /d.o.f . of 
the fit is 326/221, and for a IT model with APEC code, 
the x 2 /d.o.f. of the fit is 546/222. Contributions to x 2 are 
on the same scale for direct comparison. 
Bottom Panel: data for the bin 11' -14' and contributions 
to x 2 for a IT model with MEKAL code, the x 2 /d.o.f . of 
the fit is 729/221, and for a IT model with APEC code, 
the x 2 /d.o.f. of the fit is 886/222. Contributions to x 2 are 
on the same scale for direct comparison. 



improve the fit, \ 2 — 5.4 for 6 d.o.f) which, together with the 
enhancement in the O abundance in the inner regions points 
towards an increase in contribution by SNII, since O is basi- 
cally produced only by this kind of supernovae. 

The ratio S/Fe, and in general M87 data, suggests an agree- 
ment with delayed deton ation models (in particu lar for the 



inner bins), as stressed by Finoguenov et al (2001). If we con- 
sider the set of theoretical values for SNII and W7 SNIa models 
the behavior of these ratio would indicate an increasing con- 
tribution by SNIa going outward to the center. We recover the 
correct behavior if we choose the WDD1 y ield and we reduce 



the S SNII yield of [Nomoto et al. (1997)| by a fact or of two 



to three, as w as already indicated by ASCA data (Dupke & 



White 2000a) 
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Figure 4- MOS abundance profiles for O, Si and Fe, ob- 
tained with a 2T model using MEKAL code for the inner 
3 arcmin and with a IT model using MEKAL code for 
the outer bins. Uncertainties are at the 68% level for one 
interesting parameter (Ax 2 = 1). 
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Figure 5. Same as Fig. | but for Mg, Ar and S. 
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Figure 6. Same as Fig. | but for Ne, Ca and Ni. 
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Figure 7. MOS abundance ratio profiles for O/Fe and 
S/Fe. Also showed are the abundance ratios predicted by 
SNe models taken by Nomoto et al. 1997: the dashed line 
refers to the SNII model, the dash- dotted line to the W7 
SNIa model, the dotted line to the WDDl SNIa model and 
the three dotted- dashed line to the WDD2 SNIa model. 



5. Discussion 

The model emerging from the ASCA and BeppoSAX data for 
the explanation of abundance gradients in galaxy clusters was 
that of a homogeneous enrichment by SNII, the main source 
of a elements, maybe in the form of strong galactic winds in 
the proto-cluster phase and the central increase in the heavy 
element distribution due to an enhanced contribu tion by SNIa, 
strongly related to the presence of a cD galaxy (Fukazawa et 
al. 1998|, |5upkc fe White 2000bj |Finoguenov et al 200C| , pe 



Grandi & Molendi 2001, Makishima et al. 2001). As a text- 



book example we can consider the case of A496 observed by 
XMM-Newton ( [Tamura et al. 200l[ ). The O-Ne-Mg abundance 
is radially constant over the cluster, while the excess of heavy 
elements as Fe, Ar, Ca and Ni in the core is consistent with the 
assumption that the metal excess is solely produced by SNIa 
in the cD galaxy. The crucial ratio for the discrimination of 
the enrichment by the two types of supernovae, O/Fe, is then 
decreasing towards the center. 

The XMM results for M87/Virgo complicates this picture. 
Together with the improved accuracy in the measure of the 
heavy elements gradients, showing the increased contribution 
by SNIa, the statistically significant enhancement of a ele- 
ments O and Mg and the constant O/Fe ratio point toward 
an increase in contribution also of SNII. A mechanism which 
enhances the contribution of SNII over the homogeneous con- 
tribution all over the cluster must have been at work in the 
center of the cluster. Although there is little or no evidence of 
current star formation in the core of M87, the O excess could 
be related to a recent past episode of star formation triggered 
by the passage of the radi o jet, as we see in cD ga laxies with 
a r adio source (A1795 cD : van Breugel et al. 198J and A 2597 
cD: Koekemoer et al. 199£ ) .nearby (Cen A: Gra ham 199c) and 



distant radio galaxies (v&n Breugel et al. 1985 van Breugel & 



De y 1993 , Bicknell et al. 200C) (for a comprehensive discussion 
see [McNamara 1999[ )" 
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To put the above idea quantitatively, waiting for a true depro- 
jection of our data, we use previous ROSAT estimate of the 
depr ojected electron density in the center of the Virgo clus- 
ter (Nulsen & Bohringer 1995) to calculate the excess mass 
of oxygen. To estimate the excess abundance we fit the inner 
bins, where we see the stronger increase in the O abundance, 
with a constant obtaining an abundance of 0.32, while for the 
outer bins we obtain an abundance of 0.21 so the excess is 
0.11. Then we estimate the oxygen mass, Mo to be Mo = 
M H , where Mb = 0.82 n e ^ (R out - R m f , 
being the bounding radii in kpc of the bins. 



CXC<: .S.s 



Grevesse & Sauval 



A yo,&Zo 
Rin and R on 

Ao = 16 and yo<@ = 6.76 x 10~ 5 from 
(1998)[]Wc obtain a rough estimate of 10' Mq which, assum- 
ing that ab out 80 Mb of star forma tion are required to gener- 
ate a SNII ( [Thomas fc Fabian 199(i| ) and Nomoto SNII oxygen 
yield, requires a cumulative star formation of 4 — 5 x 10 7 Mq. 
This star formation is in agreement with a burst mode of star 
formation (~ 10 7 yr) at rates of ~ 10 — fflMay r 1 , as it is 
observed in the CD galaxies of A1795 and A2597 flMcNamare] 
199S|)~1 
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